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The electrodisintegration of the deuteron is considered within a relativistic model of nucleon-nucleon inter- 
action based on the Bethe-Salpeter approach with a separable interaction kernel. The exclusive cross section 
is calculated within the impulse approximation under various kinematic conditions. Final state interactions 
between the outgoing nucleons are taken into account. The comparison of nonrelativistic and relativistic cal- 
culations is presented. Partial-wave states of the neutron-proton pair with total angular momentum J — 0,1 
are considered. 



1. INTRODUCTION 

The electrodisintegration of the deuteron is a useful 
instrument which makes it possible to investigate the 
electromagnetic structure of the neutron-proton (np) 
system. Many approaches have been elaborated to de- 
scribe this reaction for last 40 years [1-5] . The simplest 
of them considered the electrodisintegration within a 
nonrelativistic model of the nucleon-nucleon (NN) in- 
teraction and outgoing nucleons were supposed to be 
free [1] (the plane- wave approximation - PWA). Those 
approaches were in a good agreement with experimental 
■data at low energies. However, further investigations 
have shown that the final state interaction (FSI) be- 
tween outgoing nucleons, two-body currents and other 
effects should be taken into account to obtain a rea- 
sonable agreement with existing experimental data at 
higher energies. Most of these effects have been consid- 
ered within nonrelativistic models [2,3]. In relativistic 
models, FSI effects could be calculated within quasipo- 
tential approaches including the on-mass-shell nucleon- 
nucleon T matrix [4,5]. 

One of fundamental approaches for a description of 
the np system is based on using the Bethe-Salpeter (BS) 
equation [6]. To solve this equation we use a separable 
ansatz [7] for the NN interaction kernel. In this case, 
we have to deal with a system of algebraic equations 
instead of integral ones [8]. A separable kernel model 
within the BS approach was not used to describe most 
of high-energy NN processes for a while since calculated 
expressions contained nonintcgrable singularities. The 
separable kernels proposed in [9-11] make it possible 
to avoid those difficulties. Using these kernels FSI can 
be taken into account in the electrodisintegration in a 
wide range of energy. The electrodisintegration cross 
section is calculated in the present paper under different 



kinematic conditions (see Table 3 below). The rank-six 
NN interaction potential MY6 [10] is used to describe 
scattered and bound 3 S±- 3 Di partial-wave states. The 
uncoupled partial-wave states with total angular mo- 
mentum J = 0, 1 ( 1 S'o, 1 Pi, 3 Pq, 3 Pi) are described 
by multirank separable potentials [11]. The obtained 
results are compared with nonrelativistic model calcu- 
lations [3] where NN interactions were described by the 
realistic Paris potential [12]. 

The paper is organized as follows. In Sec. 2., the 
three-differential cross section of the d(e, e'p)n reaction 
is considered in the relativistic impulse approximation. 
The used BS formalism is presented in Sec. 3. The de- 
tails of calculations are considered in Sec. 4. Then the 
obtained relativistic results are compared with nonrela- 
tivistic ones [3] and experimental data of the Sacle ex- 
periment [13, 14] in Sec. 5. 



2. CROSS SECTION 

When all particles are unpolarized the exclusive elec- 
trodisintegration of the deuteron d{e, e'p)n can be de- 
scribed by the differential cross section in the deuteron 
rest frame - laboratory system (LS) - as: 
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where CMott = (a cos |/2_E e sin |) is the Mott cross 
section, a — e 2 /(4tt) is the fine structure constant; M d 
is the mass of the deuteron; q = p e — p' e = (w, q) is 
the momentum transfer; p e = (E e ,l) and p' e = (E' e ,l) 
are initial and final electron momenta, respectively; Qf e 
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is the outgoing electron solid angle; 9 is the electron 
scattering angle. The outgoing proton is described by 
momentum p p (E p - 

nucleon) and solid angle Sl p = (6 P , (f) where 9 p is the 
zenithal angle between the p p and q momenta and (f> is 
the azimuthal angle between the (ee') and (qp) planes. 
Factor n = q 2 / s can be calculated through the np pair 
total momentum P squared: 



the relative np pair momentum in LS. To perform the 
integration, momenta p, q and the deuteron total mo- 
pjf+m 2 , m is the mass of the mcntum K (0) = (M d , 0) in LS are written in CM using 

the Lorenz-boost transformation C along the q direc- 
tion. The np pair wave function ipsM s 1S transformed 
from CM to LS by the corresponding boost operator A. 
A detailed description of ipsM s > the nth nucleon interac- 
tion vertex , the propagator of the nth nucleon , 
and the deuteron vertex function L M can be found in 
our previous works [8, 17]. 



s = P 2 = (p p +p n ) 2 = M 2 + 2M d LJ + q 2 , 



(2) 



defined by the sum of the proton p p and neutron p n 
momenta. The photon density matrix elements have 
the following form: 




/0 -tan 2 - + — 7° --9— (-1) 
'++~ tdn 2 + 2g 2 ' +-~ 2g 2 ' ( ) 

where Q 2 = —q 2 is introduced for convenience. The 
hadron density matrix elements 

W xx ,=W^ x e v x „ (4) 

where A, A' are photon helicity components [15], can be 
calculated using the photon polarization vectors e and 
Cartesian components of hadron tensor 



W v> = \ E \<np:SM s \j^\d:lM>\ 



(5) 



Sd s nSp 



where S is the spin of the np pair and Ms is its pro- 
jection; Sd, s n and s p are deuteron, neutron and proton 
momentum projections, respectively. The matrix ele- 
ment < np : SMs\jn\d : 1M > can be constructed ac- 
cording to the Mandelstam technique [16] and has the 
following form in LS: 

/J4 CM 
T^V x ( 6 ) 

Sp{A(£- 1 )^5M s (p CM ,p* CM ;P CM )A(£)rW( g ) x 

- (-1)> - f ) r« („ + (-irf ; K {0) ) } 

within the relativistic impulse approximation. The sum 
over n = 1, 2 corresponds to the interaction of the vir- 
tual photon with the proton and with the neutron in 
the deuteron, respectively. The total P CM and relative 
P*cm momen t a f t ne outgoing nucleons and the inte- 
gration momentum p CM are considered in the final np 
pair rest frame - center-of-mass system (CM), p denotes 



3. SEPARABLE KERNEL OF NN 
INTERACTION 

The outgoing np pair is described by the T matrix 
which can be found as a solution of the inhomogeneous 
Bethe-Salpeter equation [6]: 

T(p', P ;P) = V(p',p;P) (7) 
+ ^3 / d i kV(p , ,k;P)S 2 (k;P)T(k,p;P), 

where V is the NN interaction kernel, S2 is the free 
two-particle Green function: 



S^(k- P) = (±P + t-mf\^-k>- , (8) 

and p (p r ) is the relative momentum of initial (final) 
nucleons, P is the total np pair momentum. 

To solve the BS equation (7) partial- wave decompo- 
sition [18] for the T matrix: 

T a p, lS (p',P;P(o)) = ^2 t ab (p' ,\p'\;p Q ,\p\;s) x (9) 



,(2) 



JMab 



(y a M(-p')UcU ® (Ucyuip)^ 

is used. Here P( ) = (v^ 0) 1S the np pair total mo- 
mentum in CM, Uc = «7 2 7° is the charge conjugation 
matrix. Indices a, b correspond to the set 2S+1 Lj of spin 
S, orbital L and total J angular momenta, p = + defines 
a positive-energy partial- wave state, p = — corresponds 
to a negative-energy one. Greek letters {a, /?, 7, 6} in (9) 
are used to denote Dirac matrix indices. The spin-angle 
functions: 



yjM:LSp(p)U C = 



(10) 



■L \ r ,S p p r JM s-iSms 

% i pi ip 2 Lm L Sm s i mi i m2 
m L msm 1 m 2 piP2 

xY LrnL (p)uz {1 \p)u% 2 {2)T (- P ) 

are constructed from the free nucleon Dirac spinors u, v. 
It should be mentioned that only positive-energy states 
with p = + arc considered in this paper. Using a similar 
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decomposition for V, the BS equation for radial parts 
of the T matrix and kernel V is obtained: 

tab(p'o, \p'\;po, \p\;s) = v ab (p' , \p'\;p , \p\;s) (11) 

+ 00 OO 

+ 7^3 X! J dk ° J k2d \ k \ v ac{p' ,\p'\;k ,\k\;s) 



cd 



-00 



xS c d(k 0l \k\;s) tdb(ko, \k\;po, \p\;s). 



To solve the resulting equation (11) a separable ansatz 
[7] for the interaction kernel V is used: 

Vab{p'o,\p'\;Po,\p[,s) = (12) 
J2 X ii( s )9i a] (PoAp'\)9f ] {Po,\p\), 

where TV is a rank of a separable kernel, gi are model 
functions, Ay is a parameter matrix. Substituting V 
(12) in BS equation (7), we obtain the T matrix in a 
similar separable form: 



where: 



tab(p' , \p'\;po, \p\;s) = 

N 

Y, T ^ s )9 [ t\pM)9f{PoM) 



n j (s) = l/(X^(s) + h ij (s)), 



(13) 



(14) 



and 



hij(s) = (15) 

gl a] (k ,\k\)g l ; ] (k ,\k\) 



(v^/2 -E k + lef - fcg 



are auxiliary functions, Ef. = Vk 2 +m 2 . Thus, the 
problem of solving the initial integral BS equation (7) 
turns out to finding the gi functions and parameters Ay 
of separable representation (12). They can be obtained 
from a description of observables in np elastic scatter- 
ing [10,11,20-22]. 

4. FINAL STATE INTERACTION 

In our previous works [10, 17] the electrodisintcgra- 
tion of the deuteron was considered within the plane 
wave approximation (PWA) when the final nucleons 
were supposed to escape without interaction. Even 
though PWA is enough to describe the electrodisintegra- 
tion at low energies, it is important to take into account 
the final state interaction (FSI) between the outgoing 



nucleons. As it was shown in other works [3], the con- 
tribution of FSI effects increases with increasing nucleon 
energies or/and momentum transfer. Therefore, the rel- 
ativistic models of the NN interaction must be elab- 
orated and FSI should be included into calculations to 
get an adequate description of the electrodisintegration. 
The first relativistic model based on a separable kernel 
approach was Graz II [20]. However, it was impossible 
in principle to calculate FSI using it [11]. To solve this 
problem new relativistic separable kernels [10,11] were 
elaborated. We apply them to the deuteron electrodis- 
integration including FSI in this paper. 

The outgoing nucleons are described by the BS am- 
plitude which can be written as a sum of two terms: 

i; SMs (p,p*;P) = ^sl s (p,p*;P) (16) 

+i^3 5 2(p; P) J d 4 k V{p,k-P)^ S M s {k,p*;P). 

The first term 

i,f Ms {p,p*-P) = {2n)\ SMs (p-P)5(p-p*) (17) 

is related to the outgoing pair of free nucleons (PWA), 
Xsm s is a spinor function for two fermions. The second 
term in (16) corresponds to the final state interaction of 
the outgoing nucleons. It can be expressed through the 
T matrix if we use the following relation: 



J d A k V(p,k;P)^ S Ms(k,P*;P) 
J d 4 fcT(p,fc;P)4°2 /s (fc,p*;P) 



(18) 



and can be rewritten as 

41>,p*;^)= (i9) 

4ttzS 2 (p; P)T(p, p* ; P) X sm s (p* ; P) , 

here (t) means that this part of the np pair wave func- 
tion is related to the T matrix. Applying the partial- 
wave decomposition of the T matrix (9) the expression 
(19) can be written as follows: 

^f Ms {p,p*-P)=Amx (20) 

C L™SMs Y L m (P*)yaM(p)(f>a,J:LS+(PO, \p\', s), 

LmJMa 



where p* — (0,p*) with \p*\ = s/4 — m 2 is the rel- 
ative momentum of on-mass-shell nucleons in CM, p* 
denotes the azimuthal angle Q p * between the p* and q 
vectors and zenithal angle (f>. Since only positive-energy 
partial-wave states are considered here the radial part 
is: 
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t>a,J:LS+{PQ, \p\] S) = ? j= /r> ; 2 _ 2 ■ (21) 



(v/i/2 - £ p 



According to definition (10) spin-angle functions y can 
be written as a product of Dirac 7 matrices in the ma- 
trix representation [8] as: 



y a M{p) = 
1 1 



(22) 



8tt AE p {E p + m) 



matrices £ a j\/ are given in Table 2. Decomposition (20) 
is considered in detail in [19]. Using definition (16) and 



a= { 2S+1 L P j} 
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Spin-angular parts Q a M (22) for the np pair; pi = 
(E p ,p), P2 = (E p ,—p) are on-mass-shell momenta, 
-Ep = \/p 2 + rn 2 ; 7 matrices are defined as in [23]. 



substituting (17), (20) into (6), the final expression for 
hadron current < np : SMs\ju\d : 1M > can be ob- 
tained. It consists of two parts. One of them: 



<np: SM s \j„\d:lM > (n) = 

1 E {A(£- 1 )xsM S (p* CM ;P CM )A(£)r(")( (? ) x 



(23) 



S (n) l' K (0) 



(-i)V - \ I r M 



(W(-ir|;K (0) )} 



corresponds to the electrodisintegration in PWA. An- 
other one: 



<np:SM s \^\d:lM>^ 
J_ ST ST n JM J 

47T 3 



(24) 



E Qm/-(f) x 

n—1,2 LmJMjL'lm' 

CO 1 27T 



CM 
P 



J d P c M J {p CM fd\ P CM \ J dcose: 

-00 -10 

s P {A(L- l )yj L ,sM J (p CM W)r < i r ) (q) 



K, 



p+i-iy 



^ £ (p cm > cm |;Q>1;s) x 

(yi/2-E p + ie) 2 -p2 

mU + (-iT^,P+(-i)%K {0) 



corresponds to the process when FSI is taken into ac- 
count. Here gi is a radial part of the deuteron ver- 
tex function T M . The part Sp{...} has been calcu- 
lated using the algebra manipulation package MAPLE. 
The three-dimensional integration over J5q m , |p CM | and 
cos 9p M has been performed numerically using the pro- 
gramming language FORTRAN. 
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Cross section (1) depending on recoil neutron momen- 
tum |p n | calculated under kinematic conditions set I of 
the Sacle experiment [13]. The notations are follow- 
ing: MY6 (PWA) (red solid line) - relativistic calcu- 
lation in the plane-wave approximation with the MY6 
potential [10]; MY6 (FSI) (blue dashed line) - relativis- 
tic calculation including FSI effects; NR (PWAc) (violet 
dotted line) - nonrelativistic calculation [3]. 



5. RESULTS AND DISCUSSION 

The differential cross section (1) is calculated un- 
der three kinematic conditions of the Sacle experiment 
[13,14] (described in Table 3) and is present in Figs. 1-3. 
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set I [13] 


set II [13] 


[141 

L J 


E e , GcV 


0.500 


0.500 


0.560 


E' e , GeV 


0.395 


0.352 


0.360 


0, ° 


59 


44.4 


25 


v GcV/c 


mill 


0.005 


0.165 


0.294 




max 


0.350 


0.350 


0.550 


#n, ° 


min 


101.81 


172.07 


153.01 




max 


37.78 


70.23 


20.81 


a o 


min 


48.79 


44.74 


33.06 


p p , GeV/c 


min 


0.451 


0.514 


0.557 




max 


0.276 


0.403 


0.306 


° 


min 


0.622 


2.54 


13.86 




max 


51.03 


54.90 


140.28 


n o 


min 


49.41 


47.28 


46.92 




max 


99.81 


99.64 


173.35 


y/s, GeV 


1.929 


1.993 


2.057 


V^i - 2m, GeV 


0.051 


0.115 


0.176 


Q 2 , (GeV/c) 2 


0.192 


0.101 


0.038 


w, GcV 


0.105 


0.148 


0.200 


\q\, GcV/c 


0.450 


0.350 


0.279 



Kinematic conditions considered in the paper. Here all 
quantities are in LS. In addition to those which are de- 
fined in the text, they are: angle 8 qe between the beam 
and the virtual photon; neutron momentum p n and 
angle 9 n between the neutron and the virtual photon 
(p p , 6 P - the same for the proton); 8 pe (8 qe ) - the angle 
between the beam and the proton (virtual photon). 




0.01 1 ' ' ' ' ' ' ' 1 

150 200 250 300 350 



p n (MeV/c) 

The same as in Fig.l but under kinematic conditions 
set II of the Sacle experiment [13]. The nonrelativis- 
tic calculation NR (FSIc) (brown dashed-dotted line) 
includes FSI effects. 
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The same as in Figs. 1,2 but under kinematic conditions 
of the Sacle experiment [14]. Two additional results are 
presented for comparison: MY6 (PWAc) (pink dashed- 
dotted-dotted line) - relativistic PWA calculation; MY6 
(FSIc) (orange dashed-dotted line) - relativistic calcu- 
lation with FSI effects; both obtained under current 
conservation condition (25). 

The calculations have been performed within the rela- 
tivistic impulse approximation for two different cases: 
when the outgoing nucleons are supposed to be free 
(PWA) and when the final state interaction between 
the nucleons is taken into account (FSI). The partial- 
wave states of the np pair with total angular momen- 
tum J = 0, 1 have been considered. The used relativis- 
tic model consists of two parts: the separable potential 
MY6 [10] for the bound (deuteron) and scattered 3 Si- 
3 Di states and separable potentials of various ranks [11] 
- for all the other partial- wave states ( 1 5o, 1 Pi, 3 Pq, 
3 Pi). The obtained results have been compared with 
the nonrelativistic model [3]. 

In Fig.l, relativistic MY6 (PWA), MY6 (FSI) and 
nonrelativistic NR (PWAc) calculations under kine- 
matic conditions [13] (set I) are practically coincide with 
each other and go close to the experimental data. The 
notation " c" in NR (PWAc) means that the correspond- 
ing hadron current satisfies the current conservation 
condition [1]: 

= °- (25) 

So, it can be concluded that relativistic and FSI ef- 
fects do not play an important role in description of 
the deuteron electrodisintegration at low energies. 

Fig. 2 considers the cross section (1) calculated un- 
der kinematic conditions [13] (set II). Here FSI effects 
decrease the resulting cross section noticeably. A slight 
difference is seen between nonrelativistic and relativistic 
results. 
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The cross section under kinematic conditions [14] is 
presented in Fig. 3. In this case PWA and FSI calcula- 
tions differ significantly. It means that the influence of 
FSI on observables increases with increasing \p n \ . To in- 
vestigate relativistic effects, two additional calculations 
MY6 (PWAc), MY6 (FSIc) have been performed when 
the condition (25) has been taken into account. The 
relativistic one-particle current within the impulse ap- 
proximation (6) does not satisfy (25). We impose this 
condition on the hadron current (23), (24) as it was done 
in [3] for a nonrelativistic model to compare our rela- 
tivistic result with it and investigate relativistic effects. 
It has been done only under kinematic conditions [14] 
since an influence of the current conservation condition 
becomes significant at high outgoing neutron momenta 
\p n \. The comparison of MY6 (PWAc) and NR (PWAc), 
MY6 (FSIc) and NR (FSIc) demonstrates that relativis- 
tic effects become significant with increasing \p n \. FSI 
effects decrease the obtained cross section. 

As it is seen from the figures, an influence of FSI 
increases with increasing outgoing neutron momentum 
\p n \. FSI effects are significant at outgoing neutron 
momenta \p n \ > 150MeV/c even at low momentum 
transfer. It has been demonstrated that the FSI con- 
tribution always decreases the value of the cross section 
(1). From present results it can also be seen that the 
FSI contribution relatively to the PWA one is approx- 
imately the same both for nonrelativistic and relativis- 
tic models. However, other effects, like negative-energy 
partial- wave states and two-body currents, should be 
taken into account for further conclusions. It should 
be also emphasized that the relativistic calculation of 
FSI has been performed within the BS approach using 
a relativistic separable potential for the first time. Re- 
cently the considered separable potentials [10,11] have 
been modified to take into account inelastic processes 
(production of various mesons) occurring in the NN in- 
teraction at high energies [22]. However, energies of 
nucleons do not reach the inelasticity threshold under 
kinematic conditions [13]. The inelasticity effects are 
nonzero under conditions [14] but they are conjectured 
to be small in this case. 
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